BLOW UP IN SEVERAL POINTS FOR THE NONLINEAR 
SCHRODINGER EQUATION ON A BOUNDED DOMAIN 



NICOLAS GODET 



Abstract. Given p points in a bounded domain of R'*, with d = 2,3, we show the 
existence of solutions of the L'^-critical focusing nonhnear Schrodinger equation blowing 
up exactly in these points. 



1. Introduction 

We consider the L^-critical focusing nonhnear Schrodinger equation posed on a bounded 
and regular domain Q, of (with d = 2,3): 

(1.1) idtu + Au = -\u\^/'^u, {t,x) e[0,T) xQ. 

We add an initial data and the Dirichlet boundary condition : 



(1.2) 



u{t,x) = 0, {t,x) £ [0,T) X di}, 
u{0,x) = Uo{x), X G O. 



Thanks to the Sobolev embedding H'^{n) ^ L~(J]), one can show that the equation (jl.ip 
is locally well-posed in the space i7^(r2)nffQ (0) : for every initial data uq G H'^{Q)riHQ{Q), 
there exists a time T G (0, +oo] and a unique function u G C{[0,T), H'^{Q) n Hq{Q,)) 
solution of (|l.ip with initial data uq. If t( is a solution of (ll.ip . the energy and the mass 
are conserved : for every t G [0, T) 

1 d 4+2d 

E{t) := -\\Vu{t)\\l,-j^\\u{t)\\J_^=EiO), 

M{t) := ||n(t)||i2 =M(0). 
Moreover, we have the following blow up criteria 

if T < +00 then ||ti(t)||j:^2 — )■ +oo when t^T. 



Note that if d = 2, the equation (jl.ip is also well posed in Hq{Q) and even globally well 
posed if the L^-norm of the initial data is smaller than ||(5||L2(]g2) (see [H [3l [T3]). 



Blow up solutions for the equation (jl.ip has been extremely studied in and we expect 
that some results remain valid in bounded domains or more generally in the setting of flat 
geometries. Among papers concerning the study of the nonlinear Schrodinger equation on 
a domain, we can mention [21 El El \TT\ . 

In [9] , Merle shows that if O is the whole space M"^ (without restriction on d) then given 
p points in W^, there exists a solution of the focusing nonlinear Schrodinger equation with 
L^-critical nonlinearity that blows up in the p points. The aim of this paper is to show 
that this result is still true if M"^ is replaced by a bounded and regular domain of M'^ with 
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d = 2,3. We can not adapt the construction of Merle because the proof cruciaUy uses the 
dispersion estimate 

3C > 0, Vt > Oyv G Li(M^), ||e**'^i;||ioo(Kd) < ^|b||Li(M.) 

which turns out to be false if we replace M"' by a bounded domain. To prove our result, we 
use the perturbation method introduced in [10] and used in [3j to treat the case of a point 
in dimension 2. Because of lack of regularity of the nonlinearity, the case of the dimension 

3 requires a change in the choice of the weighted space where we show the property of 
contraction. 

In the following, we denote by Q the unique ([SI [H]) radially symmetric and strictly 
positive solution of 

-AQ + Q = IQI^/'^Q 
and satisfying the exponential decay at infinity : 

Va G N",3C7„ > 0,3Da > 0,Vx G R'^,\d'^Q{x)\ < Cae"-^"!^!. 

2. Statement and proof of the result 
Now we state the theorem of the paper. 



p 



Theorem 2.1. Let Q be a bounded and regular domain o/M (with d = 2,3) and xi, . . . ,x. 
p distinct points in il.. Let ipi,...,ipp G C^{il.) with disjoint supports and such that 
< v^fc ^ 1 o-i^d ifk = l near x^. Then 

(1) There exists Aq > such that for every A > Aq, there exists a time Tx > and 
ux G C{[0,Tx), H'^{0,) n Hq{0,)) such that the function hx defined by 

is a solution of Moreover, 



3 7>0,C>0,VA>Ao,VtG [0,Ta), \\ux{t)\\H2^n) < Ce ^(^a-). 

(2) For A > Ao, the solution hx blows up in at time t = Tx in the points xi, . . . ,Xp 
with speed {Tx — t)~^ . More precisely hx verifies 

(i) for R>0 small enough, for every k, \\h\{t)\\^2(B(xk,R)) WQWl^Rd)^ 

(ii) for all t G [0,rA), \\hx{t)\\L2(n) = VpWQWl^r"), 

p 

(Hi) \hx(t)\'^ — > ||Q||/,2(]Kd) (^Xfc in the sense of measures, 

k=i 

Scheme of the proof. If i7 = M"^, we know an explicit solution m of (jl.ip which blows 
up in Xi. Next we consider the function u = ipiui + . . . (fpUp where ipk is a cut-off function 
near Xk- Therefore, u is a function which vanishes on the boundary and has the same 
behavior than near j;^ because of the cut-off functions. Thus, u blows up in the points 
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xi, . . . , Xp. However, u is not a solution of (II. ip but we shall show the existence of a rest r 
such that M + r is a solution of (II. ip and keeps the behavior of u when t tends to the blow 
up time. For this, we will impose that r{t) tends to at the blow up time. To prove the 
existence of the rest r, we perform a fixed point argument in a suitable weighted space. 



Proof. For T > and A > 0, we introduce 



/ X- Xk 

[x{T-t] 



r\it,x) = ^ipk{x)rx{t,x). 



k=l 



For every k,r^ is a solution of (jl.ip on M*^. This solution blows up in H^(R'^) at time 
T and in x^. We seek a condition on ux for which the function hx := rx + ux is a solution 
of (jl.ip on i}. We have the equality 

p 

idthx + Ahx = Y, (^A A<pfc + 2Vipk ■ Vrl - ^fclr^l^/'^rt) + {idt + A)nA. 



fc=i 

Thus, hx is a solution of (jl.ip if and only if ux satisfies 

p 



{idt + A)nA 



\rx + ux\^''^{rx + ux) 



So + S{ux) 



k=l 



(rx ^Vk + 2V^fe • Vri - ^k\u\\^''^u\) , 



where we denote 



Sn 



E 

k=l 

Y^k 

k=l 



fkri 



4/d / p N 



4/(i 



\k=l 



k=l 



\k=l 



U 



+ Y'^k1 



k=l 



4/d 



U 



k=l 



To not perturb the behavior of rx at the blow up time, we impose the condition u{t) tends 
to when t tends to T ; this leads to consider the integral formulation 



(2.1) 

We introduce 



uit) 



hit) = i 

i{um = i 



3^(*-^)^(So(s) + 5(t/)(s))ds. 



At-s)^S{u{s))ds. 



To estimate the terms Iq and / we begin with the following lemma. 



4 N. GODET 

Lemma 1. There exists a constant C > such that for every u,v ^ 
- \v 

|4/d 

< C\\u 

Proof of Lemma Ul (i) ■ We use the Taylor formula 



{i) lll'ul^/'^u - |t;|^/'^t;||i2(n) < C\\u - v\\l2(^<^) {\\u\\L°-(n) + \\v\\L°-{n))'^^'^ 
[ii) \\\u\'^/'^u - \v\'^''^v\\hi(q) < C\\u - v\\H2(n) (||u||^2(t^) + ||n||^2(f^))^/'^ , 
{in) |||'u|^'''^n||^^2(n) < C||n||J^2^(5J. 



(2.2) f{u)- f{v) = {u-v) dj{tu + {l-t)v)dt + {u-v) d^f {tu + {I - t)v)dt 

Jo Jo 

with the complex function f(z) = \z\^^^z. The computation of the derivatives of / shows 
that 

\dJ{z)\ + \chf{z)\<C\z\^/''. 
We deduce from (j2.2p that for t; G C 

' < C\u-v\ {\u\ + \v\)^^'^. 



\u\'/'iu-\v\'/''v 



Then, we apply this inequality to the functions u and v, integrate and the conclusion 
follows by using Holder inequality. This prove the first point. 

(ii). First, using (i) and the Sobolev embedding H^{n) ^ L°°(0) we get 

- \v\^^'^v\\l2 < C\\u - v\\h2 {\\u\\H2(n) + \HH2in)T^'^ ■ 



(2.3) 



Next a direct computation shows that 
(2.4) di{\u 



u . 



According to the last identity, we can split di{\u\^/'^u) — di{\v\'^^'^v) into two terms. These 
two terms are treated in the same way. Let us for instance treat the second. We may 
write 

{diu)\u\^/''-^u^ - {div)\v\^/''-^v^ = A,+ A2, 

where 

Ai = ^^{u-v){\u\^/'^-^U^) , 

A2 = div{\u\^'''-^u^ -\v\''l''-^v^) . 
But by Holder inequality and again the Sobolev embedding of into L°°, 



^ II 4/d 

< ||n - 1)11^211^11^2 • 



For ^2, we first use (|2.2p to get 



<C\u-v\{\u\ + \v\fl'^-^. 



^^/d-l 



Hence 

11^2 ||l2 < \\v\\h1 \\u - v\\l^ (||^^||l°° + \\v\\l° 

< \\u-v\\h2 {\\u\\h2 + WvWn^f'^ . 
Summing the estimates on Ai,A2 and (j2.3p . we obtain the second point. 
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(iii). The norm || • \\u2 is equivalent to the norm || • \\i^2 + ||V^ • \\i2. Using again a 
Sobolev embedding, we deduce that 

(2.5) ||H^NlL^<||</.'+'<C||<f\ 
Moreover, deriving the relation (j2.4p . we get 

(2.6) + ^ ( + ^) ^^u^ju\u\'''^-^u + diju\u\''/'^~^u^ ^ ~ ^^^^i^l^l'^'^"^^^ 
Then using an Holder inequality on each term, we get 

By the embeddings L"^ and L°°, we may write 

(2.7) ||a.,(|n|^/'^n)||z.2 <C||tx||^+//'. 

Gathering (j2.5p and (j2.7p . we obtain the third point. □ 



Since the supports of are disjoints, 5'o(t) is zero near the points x^. Therefore there 
exists r > such that 

Vt G [0, T), Vx G \ U^^,;B(xfc, r), 5o(i, a:) = 0. 

Using the exponential decay of the ground state Q and its derivatives, we get the existence 
of C, D such that for every k, t, A, a with |a| < 3 : 

Therefore, Lemma [1] (iii) and the structure of algebra of H^{Q. \ U^^^i?(xfc, r)) allows us 
to write 

(2.8) \\Sm\H'^m = l|5o(t)|lH2(n\uLi5{x.,0) ^ ■ 
Now, we can introduce the following metric space 

Et =\ue L°^{[Q,T),H'^ r\Hl), sup (e^^) \\u{t)\\L2) + sup {e^^)\\u{t)\\H2) < l\ 

y 0<t<T ^ ' 0<t<T ^ ^ J 

equipped with the distance 

d{u,v) = sup (e^(^\\u{t) -v{t)\\L2) , 

0<t<T ^ ' 

where a is a real number such that < a < min(l, | — 1). We are going to perform the 
Banach fixed point argument in Eq^ to prove that the map defined by 

(2.9) ^{u)(t) = i r e*(*-^)^ (So(s) + S{u){s)) ds 



has a fixed point. For this, we show that for T small enough and A big enough, <I> sends 
Et into itself and is a contraction. 

First, we prove that {ET,d) is a complete metric space. It suffices to show that Ej- is 
closed in the complete metric space E = {u £ L°°{[0,T), L'^),exp{j^7^—j-^)\\u{t)\\i^2 < 1} 
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equipped with the distance d. Let (u„) be a sequence in Et tending to u ^ E for d. 
Since Vn{t) := exp{j^jr-j^)un{t) is bounded in L°°{[0,T), H"^ n Hq), we can extract a 

subsequence tending to f € L°°([0, T), i/^ n//g) for the weak*-topology. Then necessarily, 
by uniqueness of the hmit in P'((0, T), v{t) = exp{j^0-j^)u{t) . Using the lower 

semi-continuity of the norm in L°°([0, T), H Hq), we have 

sup fe^^(^||ii(t)||j:^2) < liminf sup (e^(^ \\un{t)\\HA , 
0<t<T ^ ' 0<t<T ^ ' 

hence taking the limit inferior in the inequality 

sup (eM2^||n„(t)||i2) + sup (e^^)\\un{t)\\u^ <\, 

0<t<T ^ ' 0<t<T ^ ' 

we obtain that u G Et- 

Boundedness. Let T £ (0, 1), A > 1 to be chosen later. We prove that $ sends E^ in 
itself. We may write 

S(u) - S(v) = \rx + v\^/'^{rx + v) - \rx + u\^/\rx + u). 

Lemma [Hand the Sobolev embedding ^ provides 

\\s{u){t) - sivmh^ < c\\u{t) - vmL^ ihmm + ii^wk^ + wrxmL-f ■ 

But using the explicit formula of rx, one can compute its derivatives to get for /c = 0, 1, 2, 



(2.10) l|VVA(t)|Uoo(n) < 

and deduce that if u,v £ Et, 



< c 



\\S{u){t) - Siv){t)\\L2 < Ce-^)d{u,v) (l + ^ 



A2(r-t)2 



By integrating this inequality, and using 

rT 



e HT-s) 



{T-s 



ds < CXe MT-t) . 



it follows 

(2.11) \\I{um - I{v)mL^ < Ce~^)d{u,v) (^r+ i 

Now we need to bound ||/(n)(i)||j|^2. The norm || • ||//2 is equivalent to || • ||l2 + HV^ ■\\h2. 
In the formula (j2.6p . there are two types of terms. The first one are those containing the 
product of two first derivatives, namely 

diudju\u\^/'^-'^u, d^udju\u\'^/'^~'^u, diudju\u\'^/'^~'^u, diadju\u\'^^'^~^u^; 

and the second are those containing a second derivative, namely 

diju\u\^/'^, d,,u\u\^/'^-^u\ 

Each term of the same type are treated in the same way. Let us begin with the first type, 
for instance the term diudju\u\'^/'^~^u'^ . We may write 

^^rK^j7%rx\^''^"^rl - di{rl + u)dj{r^ + u)\rx + uf'^-\rx + uf < Bi + B2 
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Bi 
Bo 



L2 



where 

[dir^dju + diudjvl + diudju) \rx + n|^/'^~^(rA + u)^ 
Using the estimate on the derivatives of r\ (|2.10p and the inequality 

,|4/«(-4^3 _ |^|4/d-4^3| < c\u - v\''/''-\ 

we have successively (with always u E E^)- 

Bi<C\\d,rxdjrx\u\^/''-'\\L2 

< 



L2 



c 



0(5 



|4/d~l 



< CA^e ^(^-*) . 
For B2, Holder inequality provides 

B2 < \\\rx + u\^/^~^ {^^rx^ju + diudjrx + d.udju) \\l2 
< \\rx + uf/t' {\\Vrx\\L^\\Vuh2 + \\Vu\\%) . 

But by interpolation between the spaces and H^, we remark that there exists /? G (a, 1) 
such that for every u G St, and for every t G [0, T), 

/35 



Hence 



^2 < C 



+ 1 



g A(T-t) _|_ g A(T-t) 



< CAe ^(^-') . 



Now we treat terms belonging to the second type. For instance the term diju\u\'^/'^ ^n^. 
We have 



where 



yijrx\rx 

Ci 
C2 



l^/'^-Vl - d,,{-x + u)\rx + uf'-^rx + uf^^ <C, + C2 



d.,-x(\rx\"''-\l-\rx + u\'''-\rx + uf) 
d,,u(\Tx + utl''-\rx + uf) 



L2 



The inequality 



;,|4/d-2^2 _ |y,|4/d-2^2 < _ ^1 (1^1 ^ |^|)4/d-l 
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applied with z = rx and w = rx + u allows to write 

Ci < ||V rx\\L°°\\u\\L^ (^IfaIIloc + \\u\\j^o. 

c 



< e h 1 



For the second bound, we may write 

C2 < \\'^'^u\\i^2\\rx + u\f/^ 
[ 1 

< Ce MT-t) __ _ + 1 



A2(r-t)2 

Summing Bi, B2, Ci and C2, we deduce for A > 1 and u G Et, 
(2.12) ||5(n)(i)||^2 <Ce~Mf^ fA2+ ^ 



A2(r-t)2 
and integrating we get 

(2.13) mnmWH^ < Ce-W^ (tX^ + i 



Moreover by ()2.8p . we obtain the bound 



,5 



||/olb2 < CTe MT-t). 

Finally, the latter estimate together with (j2.13|) and (j2.1ip applied with v = gives for 
A > 1 and u S Et 

(2.14) sup (e^(^\\<^{u){t)\\L2) + sup (e^<^ \\<^{u){t)\\H2) < C ( + \] . 

0<t<T ^ ' 0<t<T ^ ^ V ^/ 

Contraction. Actually, we have already proved the property of contraction of during 
the proof of the boundedness. Indeed, the estimate (|2.1ip provides for every u,v & Et, 

(2.15) d mu), <^{v)) <c(t + j^ d{u, v). 

Conclusion. First ^{Et) C L°°{[0,T), H"^ ^ H^)- Indeed, the bound (l2Ji]) gives 
^{Et) C L°^([0, T), ff2) a]2d it remains to verify that $(n)(t) G Hq almost everywhere. 
For this, it suffices to prove that for u G H"^ n H^, S{u) G Hq. But if u G n Hq, we can 
appoximate u in H^-novm by a sequence m„, G Cg". And by Lemma [T] (ii), we obtain the 
convergence in i?^-norm of |n„|^/'^u„ to |ii|^/"'n and this shows that G Hq. Hence 

S{u) G and ^{Et) C L'^{[0,T), H'^ HH^). Moreover estimates ([2Til) and (|2T5D prove 
that we can choose A big enough and T > 0, depending on A, such that ^{Et) C Et and 
for every u,v £ Et, 

(2.16) d ($(ti), <^{v)) < ^d (n, v) . 

Thus, we can apply the Banach fixed point argument with the function <I> and this 
proves the existence of the rest u G Et satisfying (j2.ip . To obtain the continuity in time 
with values in H^, we use the integral formulation satisfied by u. Estimates (|2.8p and 
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(f232]l show that the map s ^ So{s) + (5(u))(s) belongs to L°°{[0,T), H^) and since u 
verifies 

rT 

~isA 



{So{s) + {S{u)){s))ds, 



u{t) = ie''^ 

we conclude that u £ C([0, T),H'^). 

Let us show the second part of the theorem, (i). For R > small enough, we remark 
that for X G [JkB{xk, R) 



k=l 



Then 



\\hxm 



Using the fact that the norm of is only concentrated in x^ with mass ||Q||L2(Kd), and 
that ||ttA(i)||L2 converges to when t goes to Tx, we obtain (i) by a passage to the limit. 

(ii) . We may write 

p 

Whximh = hximh + hkrimh^n) + 2Rei{rx{t),ux{t))L^). 

k=l 

The first and the latter term goes to when t tends to T\ because ||itA(i)||L2 tends to 0. 
And the second term converges to p||Q|||2 thanks to the property of concentration of 
near Xi. 

(iii) . Let ip he a continuous function with compact support. Then if we denote 

r ^ 



fe=l 



we have 

m\ < c 



\rx{t,x)\'^^P{x) - Q\x)YHxk) 



k=l 



dx + C I \u\{t,x)\^dx 



< 



k=l 



rl{t,x) ij{x)-Q'^{x)^{xk] 



dx + Ce ^(^A-t) 



fc=i 



{I - ^l{x))il:{x) rlit^x) 



dx 



< 



fc=i 



\rl{t,x) ip{x)-Q^{x)i;{xk) 



P ^ 
dx + CY^ 



rl{t,x) 



dx 



+ Ce ^(^A-t). 

But converges to ||(5||^2md\'^Xfc when t — T;^ so the first term goes to 0. And the 



second one as well by the well known properties of r^. 
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(iv). We have the equahty 



k=l 



k=l 



Remarking that ||VtiA(t)||L2(Q) decays to when t goes to Tx and || Yll:=i^x(.'^)'^^k\\L^{n) 
is bounded, we get the equivalence 



l|V/lA(t)||L2(n) 

But the have disjoint supports so 

p 



p 

E 

k=l 



k=l 



and for all k, 



^k^Ait) 



L2(Q) 



L2{C) t^Tx 




^kVrlit) 



fk'^rlit) 



L2(n) 



L2(f7) 



1/2 



Vrl{t) 
1 



t^Tx X{Tx 

We obtain (iv) by summing these equivalences. 
Let us give some remarks about Theorem 12.11 



L2{n) 
^l|VQ||L2(Kd). 



□ 



Remark 1. The existence of blow up solutions for the L'^ -critical focusing nonlinear 
Schrddinger equation remains true if we replace the hounded domain of by a Rie- 
mannian manifold of dimension d = 2,3 which is locally isometric to an open subset of 
M'^ near the blow up points ; that is the case of the flat torus T*^ = R'^/Z'^. 

Remark 2. For d = 3,4, one can construct solutions of the following L"^ -supercritical 
equation posed on T'^ 



(2.17) idtu + /S.fdU = —\u\<i-^u, 

blowing up on the union the p circles. Indeed, for xi, . . . ,Xp G T'^"^, by Theorem \2.1[ there 
exists a solution u £ C([0, T), i?^(T'^~^)) blowing up in the p points. Then, we consider 
the function v E C([0, T), i?2(T'^)) defined by v{t,a,b) = u{t,a) for a ^ T'^'^ ,b & T. Then 
V is a solution of 1^2.17^ and blows up in the p circles {x^} x T. Note that the blow up on 
a sphere for the supercritical equation has been studied more precisely in M" (see |12j ]. 

Remark 3. One can also interest in the case where the equation is posed in dimension 
greater than 3. However, in this case, the nonlinearity is not regular enough to perform 
the same proof than in the case d < 3. Indeed, we need to solve the equation in a space 
included in with s > d/2 > 2 to get the embedding into L°° but if d > 4, we can not 
derive the nonlinearity more than two times. 

Acknowledgements. I would like to thank Nikolay Tzvetkov for introducing this subject, 
for his advices and his helpful remarks on my work. 
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